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ABSTRACT 


The concept of gravity tectonics is applied to reveal the major clue as to the 
conditions which result in the correspondence of seismic and tectonic gaps in the 
mantle. An asymptotic theory is developed for the calculation of the thrust and 
moment when a descending lithospheric plate encounters resistance to its down- 
ward motion in the mesosphere. Dynamic analysis falls into two parts: (1) de- 
riving equations for forces in the descending lithosphere, (2) deducing moment 
distribution which causes the detachment of lithosphere. For the analysis of 
forces a mathematical theory of shells is given. In order to determine the de- 
tachment mechanism, solutions of equations are obtained by asymptotic integra- 
tion. It is found that a thrust coupled with a moment due to gravitational 
forces generated by density contrast may play a key role in the initial detach- 
ment of a piece of descending lithosphere. The results are in agreement with 
the observed seismic gaps beneath South America, Toga-Fiji, New Zealand and 


New Hebrides regions. 
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INTRODUCTION 

Observations of the variations in attenuation and velocity of seismic 
waves beneath island arcs (Barazangi, et al. , 1973; Isacks and Barazangi, 1973; 
Pascal, et al., 1973; Stauder, 1973) have revealed remarkable gaps in seismicity 
These gaps are of interest because of the implication that portions of the litho- 
sphere can break off from the descending plate and exist as isolated slabs in the 
mantle. Although the existence of isolated lithospheric slabs is a natural expla- 
nation of much of the data on which seismic gaps are based, the dynamic aspects 
of detachment of these slabs are unknown. The present paper is devoted to study 
ing the tectonic setting of seismic gaps beneath island arcs. 

In this paper a model of gravity tectonics (Hales, A. L. , 1969; Jacoby, 
1972) is developed to explain conditions under which the correspondence of gaps 
in seismiciiy and tectonics may occur. Computation of the gravity effect of the 
downgoing lithosphere beneath island arcs (Elsasser, 1969; Isacks and Molnar, 
1971) is complicated by its material properties. In order to obtain the first 
approximation, the density contrast between the descending lithosphere and 
surrounding asthenosphere is assumed to be uniform (Oxburgh and Turcotte, 
1970) and the underthrusting plate is assumed to be elastic to maintain its shape 
(Sykes, 1966; McKenzie, 1969; Liu, 1973; Watts and Talwani, 1974). Analysis 
of the edge effect of thrust and moment (Love, 1944; Girkmann, 1956) due to 
gravitational body forces reveals a major clue for the initial detachment of de- 


scending lithosphere. 
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EQUATIONS OF THE PROBLEM 

Several studies on the thermal regimes of downgoing slabs have been 
carried out using different models and techniques. (Turcotte and Oxburgh, 1968, 
McKenzie, 1969, 1970; Minear and Toksoz, 1970; Hasebe, et al. , 1970; Griggs, 
1972; Toksoz, et al. , 1971, 1973) In addition to heat conduction from the sur- 
rounding mantle, the slab is heated by internal heat sources. These consist 
of radioacticity, phase changes, shear heating and adiabatic compression. 
Ideally, we would like to be able to specify the properties of the mantle 
and solve the time -dependent equations of deformation for the descending slabs, 
including the effects of temperature and pressure dependent parameters 
and internal heat sources. This task is analytically impossible and presents 
formidable problems numerically. Furthermore, it may well be that unless 
simpler models are first understood, complex cases would not be interpret- 
able in terms of the contribution to the dynamics that each case provides. It 
seems necessary, and perhaps even desirable, to consider relatively simple 
model problems. 

Assumptions and simplifications, common to all model problems, are 
required to obtain a tractable set of governing equations of detachment for the 
descending slabs in the mantle. We seek to generate models sufficiently simple 
to be solvable but which retain the main elements of the dynamics of the litho- 
spheric detachment. One of the serious assumptions in tectonic plate theory is 
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the rigidity of the plate. Its validity needs to be justified by the first order of 
approximation of the dynamical principle of deformable bodies. In this model, 
the descending slab is considered to be elastic, but the role of internal heat 
sources and phase boundaries are not considered. The temperature and pressure 
influences on the mechanical properties of the slab are not included within the 
scope of this paper. Linearized density in the slab and in the asthenosphere is 
assumed. The mesosphere, lying below 600 km depth, is assumed to be more 
dense than the asthenosphere, but no specific values for density are required. 

The assumed density contrast of about 0.05 gem* 3 between the slab and its sur- 
rounding asthenosphere is the dominant feature of this model. These assumptions 
and simplifications, which impose limitations in representing the geophysical 
problem, lead to analytical solutions of the model that can be confirmed by ob- 
servations of the tectonic plates . 

Let us consider an element cut from the outer shell of the earth by two 
adjacent meridian planes and two sections perpendicular to the meridians (Fig. 1). 
In a axisymmetric shell the following resultant forces and moments per unit length 
occur: (1) in the cross-section <P = constant: N^, and (2) in the cross- 
section d = constant: N 0 and . The stresses can be reduced to the resultant 
force N^r L d^» and resultant moment M^rjd0. The side of the element perpendic- 
ular to the meridians which is defined by the angle <j> is acted upon by normal 
stresses which result in the force N^r 2 Sin0d0 and the momentM^r 2 Sin0d0 and 
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by shearing stresses which reduce the force Q^r 2 Sin0d0 normal to the shell. 
The three equations of equilibrium are (Love, 1944; Girkmann, 1956) 


— (V 0 )-rV, cos <p - Q(4,r 0 - 0 


Vo +N 0ri sin^+ — «Vo> = 0 


(1) 


“ 0V o ) ' M 0 r i cos 0 " <Vo r i " 0 

in which the resultant forces , N 0 , shearing force Q 0 and resultant moments 
M 0 and Mg are five unknown quantities. The number of unknowns can be reduced 
to three if we express N„, N 0 , M 0 and Mg in terms of the components V and W 
of displacement (Fig. 2). The strain components of the middle surface of the 
element are 


1 dV W 

i'j d<j> t'j 

from which, by Hook's law, we obtain 

Eli 


V W 

€ e ~ T cot 0 " ~ 


( 2 ) 


% = 


1-^2 

Eh 


Nfl = T 

y 1 - r;2 


i. (“X _ w) (V cot 4 > - W) 
r, \d0 / r 2 


! * /dV 

— (V cot 0 - W) h I — — -w 


r i \ d 0 


(3) 


Where E is modulus of elasticity, h is the thickness of the shell and v is 


Poisson's ratio. 


To obtain similar expressions for and M 0 , we consider the changes 
of curvature of the shell element. Considering the upper and lower side of that 
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element, the initial angle between these two sides is dp. Because of the dis- 
placement V along the meridian, the upper side of the element rotates with 
respect to the perpendicular to the meridian plane by the amount V/r r As a 
result of the displacement W, the same side further rotates about the same axis 
by the amount dW/(r, d p). Hence, the total rotation of the upper side of the ele- 
ment is V/r, +dW/(r 1 d<£). For the lower side of the element the rotation is 

V dW d / V dW 

-f b I — + 

r, r, dp dp\r l d p 

Therefore, the change of curvature of the meridian is 



1 d /V 

* " r, dJb, 


+ 



(4) 


To find the change of curvature in the plane perpendicular to the merid- 
ian, we observe that the normal to the right lateral side of the element makes 
an angle tt/2-cos p d0 with the tangent to the y-axis. Therefore, the rotation of the 
right side in its own plane has a component with respect to the y-axis equal to 


V dW 

- + co$ A 

x 1 r, d p 

This results in a change of curvature 

V _ dW\ cot p 

‘1 r . d V r 2 
Using equations (4) and (5), we have 

= -DiX 0 + 


x#= (t + 


(5) 


= -D 


1 d / V dW\ v /V dW\ 

— — f — + ] + — ■ [ — + — — | cot p 

r, dp \r t fj dpj r 2 ^r, x^dpj 


( 6 ) 
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Mg D(Ag + v\q) 

= -D 


_V dW \ cot nf> j> £ / V J^L\ 

r i r i d ^/ r 2 r i cI 0\ r i r t d0/ 


(6 Cont’d.) 


where the flexural rigidity is 


D = 


Eh 3 

12(1 -j/2) 


Substituting equations (3) and (6) into equation (1), we obtainthree equa- 
tions with three unknown quantities V, W and ♦ Simplification of these equa- 
tions can be obtained by transformation to new variables. 


TRANSFORMATION OF EQUATIONS 

By using the third equation in (1) the shearing force can be elimi- 
nated and the three equations reduced to two equations with the unknowns V and 
W. Considerable simplification of the equations can be obtained by introducing 
new variables. As the first of the new variables, we take the angle of rotation of 


a tangent to a meridian. We define this angle by 


. V dW 

<f> = — + . 


r, d 0 


(7) 


As the second variable we take the quantity 

u = r 2 (8) 

To simplify the transformation of equations to the new variables, we may 
consider the forces in the portion of the shell above the parallel circle defined by 
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the angle ip. They are governed by 

r o sin 0 dd + r Q cos 0 d 8 - 0 

from which 


. , cot 0 

N * = u 
r 2 


Substituting equation (9) into the second of equations in (1), we find 

M 1 du 

N 0 = — 


( 9 ) 


rj dtp 


( 10 ) 


Thus, and are both expressed in terms of U which is dependent on as 
defined by equation (8). 


To establish the first equation connecting and U, we use equation (3) 
from which we obtain 

dV 




Vcot^-W = (N r .N^) 


Eh 


(ID 


( 12 ) 


By eliminating W from equations (11) and (12), the result is 

dV 
d <p 


dV 1 

-V cot<p = ~ [(r t + v r 2 )N^ - (r 2 + v r, )N 0 j 


( 13 ) 


Differentiation of equation (13) with respect to <p gives 
dV V dW d 


cot (p 


d <p sin 2 cp dip dip 


Et (N »- N d 


( 14 ) 
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By eliminating dV/d0 from equations (13) and (14), the result is 



dW 

d0 


cot 0 
Eh r, 


Krj + t 'r 2 )N 0 -(r 2 +^r 1 )N 6 ] 


1 

Eh fj 


d 

d0 


[r 2 (Ng - v N^)] 


(15) 


Substituting equations (9) and (10) in (15), we obtain the following equa- 


tion relating to <i> and U. 


r 2 d 2 U 1 

r. 2 d0 2 r. 


— {— I + — cot 0 

d<A r i/ r l 


r 2 dh 


rj h d0 


dU 

d <p 


(16) 


I /r. v dh \ 

1 -i cot 2 p + cot p ~ v ) U = Eh 'f> 


r i \ r 2 


h dp 


The second equation relating to <t> and U is obtained by substituting 
equations (6), (7) and (8) in the third of the equations in (1), In this way we find 


r 2 d 2 U 


dp' 


+ — 


d /rA r 2 r, dh 

— + — cot 0 + 3 — 

rjh d0 


d <H r i/ r, 


dU 

d0 


(17) 


1 / , 3v dh \ U 

/ cot 2 0 cot 0 + V j <t> 

r, \r 2 h dp j D 

Therefore, the problem of the membrane tectonics (Turcotte, 1974) of the under- 
thrusting shell is reduced to the integration of equations (16) and (17). For the 
case of constant thickness, the terms containing dh/d0 as a factor vanish, and 
the derivatives- of the unknowns <!> and U in both equations have the same coeffi- 
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cients. By introducing the notation 


f r 2 d2 

£(X) — + 

\r, 2 d0 2 


1 P d /r,\ r, d 

Tj 1 d0 \r,y r L d 0 


GOt 2 0 


equations (16) and (17) can be represented in the following simplified forms: 


£ (U1 + — U = Eh 

r i 


v U 

£($) <I> - -- 

r, D 


Performing the operator £ on equation (19) gives 


££(U) + v £P-j = Eh £ (<f>) 


Substituting equation (19) into (20), 


v U v r v ”1 U 

m =— 4>-- = -tt- U(u>+- u -- 

Tj D ij Eh I r ( J D 


we obtain 


/U\ v v 2 Eh 

££(U) + r£[ — \ - — X(U) r U = - — U 


, r i/ r i 


If the radius of curvature r 1 is constant, equation (22) reduces to 


/Eh v 2 \ 

“ (u)+ (t— r=° 
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ASYMPTOTIC EXPRESSION OF Q0 

By introducing, instead of the she aring force Q0> the new vnrinble 

T - sin' /2 0 (24) 


equation (23) becomes 

d 4 r d 2 r 

+C, r 

c!0 4 *" 00" 


dr 

+ C, + (4 A 4 + C ) t 

d<t> 


0 


(25) 


in which 

63 9 9 

C + + — 

° 1 6 sin 4 0 8 sin 2 0 1 0 

3 cos 0 

C) = __ 

sin ' 1 0 

3 5 

C 2 = r- + — 

2 sin 2 0 2 

! - v 1 /l2a 2 

X 4 (— 7 - + 1 

4 V h 2 

In the derivation of equation (25), r { = a is assumed. For the case of under- 
thrusting plate, a s 10 3 km and h ^ 10 2 km, a/h can be regarded as large. 
Therefore, the value of 4 \ 4 is very large in comparison with the coefficients 
C 0 , C { and C 2 , provided the angle 0 is not small. Since we shall be interested 
in moments in the leading part of the underthrusting plates in the mesosphere 
where 0 =s 7 and 7 is not small, we can neglect the terms with the coefficients 
C 0 , Cj and C 2 . In this way we obtain the equation 

d 4 r 

+ 4 X 4 7 = 0 

d0 4 


( 26 ) 
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(27) 


The general solution of equation (26) together with equation (24) gives 
= sin 1/z 0 [e X< ^ (k, cos X0 + k 2 sin X0) 

+ e~ x< ^ (k 3 cos X0 + k 4 sin X0)] 
where kj , k 2 , k 3 and k 4 are the constants of integration. They must be deter- 
mined from the conditions at the leading part of the descending plate. Since the 
moments produced by forces on the edge of a shell decrease as the distance from 
the edge increases (Love, 1944; Girkmann, 1956), it is permissible to take only 
the first two terms in equation (27) and assume 


= e x< ^ sin' /2 0 (k s cos X0 + k 2 sin X0) 


(28) 


Similar mathematical analysis shows that 4>has the same oscillatory 
character. 


APPROXIMATE SOLUTIONS 


As a basis of an approximate investigation of the bending of the under- 
thrusting plates, we take equations (19) and (20). For this purpose, these equa- 
tions can be written as follows: 


d 2 Q 0 d Q0 

+ cot 0 (ooU 0 - v) Q, = Eh 

d0 2 d0 9 

d 2 4> d<h „ a 2 Q rh 

— + cot 0 (cot 2 0 + v) = - 

d0 2 d0 D 


(29) 

(30) 


and have the same oscillatory character as shown in equation (28) 
and are damped out as the distance from the leading edge of the plate increases. 
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Since X is large, the derivative of equation (28) is large in comparison with the 
function itself and the second derivative is large in comparison with the first. 

This indicates that a satisfactory approximation can be obtained by neglecting the 
terms containing the functions % and <t> and their first derivatives in the left-hand 
side of equations (29) and (30). Therefore, they can be replaced by the following 
simplified system of equations 



= Eh 

(31) 

dd 2 



d 2 <f> 

a 2 

(32) 

II 

1 ^ 

1 ^ 



By eliminating ^from these equations, the result is 

^+4^0, = 0 (33) 

dd 4 ^ 

where 

a 2 

(3 4 = 3(1 -v 2 ) ~ 

Ir 

The solution of equation (33) for our problem is 

= K, cos ($4> + K 2 sin j 3<p (34) 

K , and K 2 are to be determined from conditions at the leading part of 
the plate. In discussing the edge conditions, it is advantageous to introduce the 
angle ^ = 7 - 0 . Substituting y-ifr for <j> in equation (34) and using the new constants 
K and a, we can represent equation (34) in the form 

= K e'^ sin (/?i j/ + or) 


(35) 
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Now, employing equation (9), we find 

= - K cot (7 - \p) sin (j 8 \|/ + a) (36) 


Substituting equation (35) in equation (31), we obtain the angle of rotation 

2j3 2 


<t> = Ke' M cos (J3^ + a) 

Eh 


(37) 


The moment can be determined by introducing equations (7) and (37) 
in (6). Neglecting the terms containing <1> in these equations for r, = r 2 = a, we 
find 




D d$ 

a dp 


K e'W sin ( ftp +a + — \ 
V2P \ 4 / 


(38) 


With the aid of equations (36), (37) and (38), the model of gravity tec- 
tonics for the detachment of descending lithosphere can readily be treated. 


CALCULATIONS 

In order to calculate and in the underthrusting plates, we must 
determine the constants K and a in equations (36) and (38). This can be done by 
applications of force analysis (Fig. 3). According to the suggestion made by 
Elsasser (1969) and Isacks and Molnar (1971), the lithospheric plate is sinking 
under island arcs and exerting a downward gravitational pull. When the sinking 
lithosphere reaches the more dense, stronger methosphere region below the 
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astheno sphere, downdip compression is evidenced in the earthquake patterns 
(Isacks and Molnar, 1971). This suggestion requires a positive density contrast 
Apof the sinking plate which is about 0.02 to 0.06 g cm ' 3 as estimated by Jacoby 
(1970) and by Oxburgh and Turcotte (1970). Such an estimation is supported by 
the observed gravity anomalies (Hatherton, 1969), Satellite gravity (Kaula, 1972) 
has the resolution to show the positive anomalies related to the oceanic trenches. 
To estimate the buoyant force on the plate, we consider the case of a plate with 
edge thrusting in the mesosphere and subject to the action of gravity. The 
vertical force is W = AApg = (JZh/Sin' 7 ) Apg , where A is the volume per unit length 
along the strike, £ is the depth of the leading edge of the descending lithosphere, 
h is its thickness and 7 is the dip of the plate measured from horizontal. The 
force F represents a push from the ridge due to the elevation of the ridge associ- 
ated with ascending convection beneath the ridge and a horizontal traction on the 
base of the lithosphere due to convection in the upper mantle. If the descend- 
ing plate is uncoupled from the horizontal lithosphere by vertical faults, no ten- 
sile force could be exerted directly on the horizontal plate from the descending 
one. If the tension is compensated by compression from the ridge (Jacoby, 1970), 
the component S = WSin 7 parallel to the inclined plate may be in equilibrium with 
a vertical supporting force G and a horizontal thrust H in the mesosphere. This 
estimate is not yet complete. When the plate continues to dip into the more vis- 
cous and denser mesosphere due to force F, the resistance is likely to grow 



rapidly, and the leading edge of the plate encounters a resistance (N^) 0 = 7 
which may be greater than S. Hence, for this case, we have two boundary 
conditions 

(N^ = 7 >RhApg 

(39) 

(M 0 ) 0 = 7 = 0 


By substituting 0 = 0 in equation (38) it can be seen that the second 
boundary condition in equation (39) is satisfied by taking the constant a — - tt/ 4. 
To satisfy the first boundary condition, we use equation (36) which gives 

ChA pg 


K> 


sin (tt/ 4) cot y 


Substituting the values of the constant K and a- in equation (38), the 

result is 

, .hApg e . 0 , sM 

v j3 cot 7 

For £ = 600km, h = 1 00 km, Ap = 0.05 gem' 3 , g = 10 3 cm sec' 3 , a = 1200 km, and 
7 = 45°, |3 = 4.45, we obtain 


(40) 


(N ^) 0 = 7 >3.5x 10’ & dvn cm ' 1 
= > 8.4 x 1 0 2 3 e"^ 3, ' ^ sin (7 - 4 > ) dyn - cm cm ' 1 


(41) 


(42) 


The corresponding values of from equation (42) are shown in Fig (4), 
In Fig (4) the maximum moment of body forces due to gravity generated by den- 
sity contrast occurs at 0 = 36°. This implies that the thrust (N^)^ _ 7 coupled 
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with the bending moment may provide a mechanism for the initial detachment 
of the descending lithosphere plate at a depth of about 460 km. Under the influ- 
ence of the bending moment, M 0 curve represents the curvature of the plate. 

This analytical result reveals that the lithosphere breaks up but does not sink 
vertically along en echelon vertical fault as shown in Figure 4. The earthquake 
distribution beneath the New Hebrides arc as obtained by Dubois (1971) is illus- 
trated in Figure 5. The horizontal extent of the deep earthquakes in Figure 5 
suggest that the detachment of the slab is controlled by the gravitational bending 
moment due to density contrast. Therefore, the gravitational bending moment, 

M may provide an explanation for the similarity among the seismic and tectonic 
gaps beneath the South American, Tonga- Fiji, New Zealand and New Hebrides 
Island arcs as observed by Isacks and Barazangi (1973), Stauder (1973), Pascal, 
et al. (1973), Barazangi, et al. (1973) and Isacks and Molnar (1971). If this 
similarity is not fortuitous, the analytical solutions of this model are probably 
relevant to the geophysical problem. Previous studies (Toksoz, et al. , 1973; 
Liu, 1973) lack the understanding of the detachment mechanism. This geophysi- 
cal phenomenon is, however, not unexpected if the elastic model of gravity tec- 
tonics developed in this paper is accepted. 

CONCLUSION 

Seismic observations provide important information of the deep-seated 
gaps along converging plate boundaries. These gaps must be examined with 
studies of the tectonic plate deformation to determine whether the results of 
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seismic observations are consistent with simple tectonic plate models. In this 
paper a mathematical model of gravity tectonics to explain and predict the seis- 
mic gaps in the mantle is developed. The results of the analysis are formulated 
by equation (36), (37) and (38). Numerical calculations lead to equation (42) and 
Figure 4. It is shown that the thrust coupled with a moment of body forces due 
to gravity on the leading portion of the underthrusting plate may play a key role 
in the initial detachment of a piece of descending lithosphere. Furthermore, the 
tectonic setting of seismic gaps presented in this paper reveals the clue as to 
the conditions under which they occur and the regional motions, forces and mo- 
ments which must be responsible for their occurrence. 
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FIGURE CAPTIONS 
Figure 1. Forces in a Plate Element 
Figure 2, Displacements in Plate 
Figure 3. Force Diagram 

Figure 4. Distribution of Moment Due to Gravitational Forces 

Figure 5. Vertical Cross Section Perpendicular to the New Hebrides Arc 

Showing Earthquake Distribution Beneath the Arc and in the Detached 
Slab 
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